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^ 1 Abstract 

Induced soft gluon bremsstrahlung associated with multiple collisions is cai- 
rn 

culated via perturbative QCD. We derive the non-abelian analog of the Landau- 
Pomeranchuk effect that suppresses induced soft radiation with formation times 
q ■ exceeding the mean free path, A. The dependence of the suppression effect on 

the SU (N) representation of the jet parton as well as the kinematic variables is 



expressed through a radiation formation factor. Unlike in QED, the finite con- 
tribution from the small x regime in QCD leads to an approximately constant 
radiative energy loss per unit length, dE/dz oc u 2 , in the high energy limit that 
is sensitive to the infrared screening scale, u, of the medium. As a function of 
the dimensionless parameter £ = Xu 2 / 'E, we show furthermore how the energy 
dependence of dE/dz evolves from the above constant for £ <C 1 to the more 
familiar (Bethe-Heitler) linear dependence for 
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1 Introduction 



Radiative energy loss of ultra-relativistic particles passing through dense matter is 
of interest not only because of its many practical applications, but also because it 
illustrates a characteristic destructive interference phenomenon caused by the finite 
formation timejl], |J, 

r(k) ~ h/AE(k) ~ 2u/k 2 L ~ 2/u6 2 , (1) 

of quanta with large four momentum, fc M = (u, k z , kj_), emitted at small angles, 
9 kj_/u, relative to the incident particle. In effect, r(k), is the minimal time 
needed to resolve the transverse wavepacket of the quanta with, Ax± ~ h/k±, from 
the wavepacket of its high energy (Eq uj) parent. Destructive interference between 
radiation amplitudes associated with multiple collisions can be expected when the 
mean free path, A, is short compared to the formation time. When r(/c) ^> A, the 
emitted quanta cannot resolve different elastic scattering centers, and the assumption 
of independent contributions from each separate scattering in the medium breaks 
down. This effect, first studied in QED and then in other field theories, is often 
referred to as the Landau-Pomeronchuck-Migdal (LPM) effect Jl], j2]. 

Interest in analogous destructive interference phenomena in non-abelian theories 
is connected with attempts to understand the weak nuclear dependence of hard QCD 
processes and the use of those dependencies as a probe of the space-time develop- 
ment of hadronization. A high energy quark or gluon passing through dense (QCD) 
matter of course suffers multiple interactions. However, in the Q — > oo limit, QCD 
factorization theorems apply that show that, as in the QED case, the soft radiation 
is emitted only from the external legs||. There are a large number of phenomena 
such as the very weak nuclear dependence of Drell-Yan yields |G| and the apparent 
nuclear independence of very high energy quark fragment at ion 0, that confirms 
this basic factorization feature of asymptotic pQCD. As shown in Ref. ||, however, 
the assumptions leading to the factorization theorem break down for sub-asymptotic 
conditions. In the Drell-Yan process, for example, final state interaction corrections 
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become important for sufficiently large nuclei with, A 2 ' 3 Aq CD > Mp Y . In deep in- 
elastic hadroporduction, nuclear dependence always exits in the xp < region and 
disappears in the xf > region only when kR/v — > 0, where k ~ 1 GeV/fm[[7]. 

In order to understand such nuclear effects more quantitatively, it is necessary to 
study induced radiation patterns associated with multiple collisions in QCD. Rough 
estimates for the magnitude of the suppression of non-abelian induced radiation have 
been made in || O, |T2L [T3, [TIJ. The estimates for dE/dz vary, however, widely 



ranging from energy independent [T(|, logarithmic energy dependent || to dE/dz oc E 



12] , [L3|| . However, a detailed study of the LPM effect in QCD has not been performed 
to our knowledge in the context of multiple collision theory, taking into account 
essential non-abelian features of the problem. The aim of this paper is to initiate 
such a study. 

While pQCD can serve only as a qualitative guide because the effective coupling, 
a s , is small only in extremely dense matter, e.g., a quark-gluon plasma at tempera- 
tures T ^> T c ~ Aqcd, h is instructive to explore its consequences in situations where 
approximations to multiple scattering can be used to simplify the problem. Additional 
motivation for this work is to compare radiation patterns due to multiple collision 



physics in pQCD with those suggested by phenomenological string modelsP, IT1, ITS 



for high energy eA,pA, and AA reactions. We show, for example, that some features 
of the induced soft gluon rapidity distributions in pQCD are similar to the hadron 
distributions predicted by multi-string Lund type models utilizing string breaking and 
flip mechanisms In particular the effective string tension, k ~ 1 GeV/fm in those 
models is analogous to constant radiative energy loss due to induced soft radiation in 
the x<Cl regime in non-abelian multiple collision theory. Furthermore, additivity of 
radiation from multiple scatterings is limited to a domain x < kX/Eq, that shrinks 
as Eq increases, just as in string flip models for multiple interactions ||. However, in 
pQCD, unlike in string models, k = dE/dz oc fi 2 ~ g 2 T 2 , is found to be sensitive to 
the infrared screening scale, /i, in the medium. Recall that other so called "string 
effects" were also found to arise naturally from interfering pQCD amplitudes for three 
jet events in e + e~ B. Finally, we note the importance of clarifying radiative energy 



loss mechanisms and interference phenomena associated with formation zone physics 
for the further development of QCD transport theory and parton cascade models |]IE| . 

Before turning to the non-abelian problem, we recall briefly how the destructive 
interference for induced radiation occurs in QED. The Fourier transform of the current 
of a classical particle undergoing m collisions at space-time points xf at which the 
four momenta change from p1_ x to pf is 

m m / u u \ 

The resulting soft bremsstrahlung spectrum is given by 
d n 

= E |eJ(fc)| 2 /(2(27r) 3 ) = -| | 2 /(2(27r) 3 ) . (3) 

The sum contains m diagonal terms, where the phase factors drop out, and m{m — 1) 
off diagonal terms involving phase factors, exp(ik(xi — Xj)). Two extreme limits are 
obvious. One is the incoherent limit where k(xi — Xj) ^> 1. In this case, the off 
diagonal phase factors tend to average to zero. This corresponds to the usual Bethe- 
Heitler limit, in which the radiated energy loss, dE/dz = —E/L r , grows linearly with 
energy with L r being the radiation length (the ratio of the mean free path to the 
average fraction of the energy radiated per collision). The other extreme limit, is the 
one corresponding to k(x{ — Xj) <C 1 for all In this case, all the phase factors 
are approximately unity, and there is an exact cancellation between adjacent terms 
in eq. (0). Only the radiation from initial and final lines contributes. This is the 
so-called the "Factorization limit" since the amplitude for soft radiation factors into 
an amplitude for multiple collisions times a current element depending only on the 
momenta of the external lines entering and leaving the reaction. Because the radiation 
intensity in the high Q 2 = — (p/ — Pi) 2 limit increases only as log(Q 2 ), the radiated 
energy loss grows only logarithmically with the number, m, of elastic collisions in the 
medium. Therefore, for a random walk leading to Q 2 oc m, the induced radiation 
energy loss per collision, dE/dm oc log(m)/m, becomes negligible for large m. In 
the general case, between these extreme limits, there is a partial contribution from 
intermediate current elements. 



The essential parameter controlling this interference effect is the ratio of the for- 
mation time (|1]) to the distance between multiple interactions. For an ultra-relativistic 
particle propagating in a straight line 

k(%i - Xj) = (u - k z )L i:j = Lij/r(k) , (4) 

where La = the longitudinal distance (time) between scattering at % and 

j. For finite deflection angles, an additional phase, — ■ r^, appears that depends 
on the transverse separation of the scattering events. Interference between current 
elements Jj and Jj occurs only if k(xi — Xj) <C 1. This requires r(k) 3> Ljj and 
k± <C Tij. The interference pattern depends also on the current correlation function, 
(Ji(k)Jj(—k)). The LPM effect in QED often refers to the specific destructive inter- 
ference pattern calculated by Migdal using the Fokker-Planck transport equation to 
solve for the probability distribution, W(x, p, p'), of scattering points and initial and 
final momenta. Monte Carlo methods have also been developed []T7 to calculate the 



development of very high energy cosmic ray air showers. The interference effect found 
in the limit that the scattering medium is much thicker than the radiation length is 
that the familiar soft 1/u bremsstrahlung frequency spectrum is transformed into 
a 1/y/uj form for uj < Eq/E lpm . Remarkably, because the characteristic energy, 
Elpm ~ 3 TeV ~ 5 ergs, turns out to be so large, this interesting prediction has yet 
to be verified quantitatively experimentally [|T8| . 

In the following sections, we calculate the induced non-abelian radiation for a 
high energy parton passing through the random color field produced by a color neu- 
tral ensemble of static partons. This idealized system is chosen to minimize the 
complications of multiple scattering theory while illustrating the essential features of 
of the non-abelian LPM effect. In section 2, we first calculate the elastic multiple col- 
lision amplitude for a spinless high energy jet parton belonging to an arbitrary SU(N) 
representation. We show how classical multiple collision cascade theory emerges from 
pQCD in both the high and low momentum transfer regions in the limit where the 
mean free path is large compared to the range of the Debye screened potentials. In 
section 3, we calculate the soft gluon radiation amplitudes in the restricted kinematic 



region range k± <C fi. This restriction limits the applicability of the results to "thin" 
plasmas, with the number of mean free paths m not large, or to very small y/rri 
gluons. The extension of the results to thick plasmas, that would be necessary to make 
contact with the Migdal (R = oo,Eq — > oo) limit 0, is not considered in this paper. 
The soft eikonal approximations used here are aimed, on the other hand, to the study 
of the breakdown of factorization in the opposite limit (E = oo, R — > oo). For high 
energy reactions involving finite nuclei, this is in fact the only physically relevant limit 
for applications of pQCD in any case. Destructive interference between the radiation 
amplitudes from jet and exchanged gluon lines is shown to limit the transverse mo- 
mentum distributions, and the non-abelian generalization of @ is derived. In section 
4, we introduce and calculate the "radiation formation factor" controlling the magni- 
tude of the suppression of induced radiation in pQCD. We emphasize the role of color 
algebra on the destructive interference pattern. In section 5, the radiative energy loss, 
dE/dz, due to soft induced radiation is estimated. Finally, a discussion of remaining 
open problems is presented in section 6. 

2 Multiple Elastic Scattering in a Color Neutral 
Ensemble 

2.1 The Model Potential 

Consider the sequential elastic scattering of a high energy (jet) parton in the random 
color field produced by an ensemble of m static partons located at Xj = (zi, x_u) such 
that z i+1 > Zi and (z i+ i — Zi) ^> where \i is the color screening mass in the 
medium. As a simplified model of multiple scattering in a color neutral quark-gluon 
plasma, we assume a static Debye screened potential for each target parton: 

^ a (q) = ^ro c , C 'e ( - iq ' Xi) /(q 2 + ^ 2 ) , (5) 

where Tf is a ^-dimensional generator of SU(N) corresponding to the representation 
of the target parton at Xj. The initial and final color indices, c, c', refer to the 
target parton are averaged and summed over when computing the ensemble averaged 



5 



cross sections. With V { a oc the ensemble averaged potential vanishes everywhere, 
{V?) oc TrT? = 0. However, since 

TrT^T) = SatSijidi/dA)^ ■ (6) 

the diagonal mean square fluctuations and the cross sections are finite. Recall that 
for SU(N) the second order Casimir, C 2 i = (N 2 — 1)/2N = Cf for quarks in the 
fundamental (dj = N) representation, while Cn = N = Ca for gluons in the adjoint 
(di = N 2 — 1 = d A ) representation. 

In this potential, each scattering leads on the average to only a relatively small 
momentum transfer gf = q zi , q_u) with each component being much less than the 
incident energy, E Q . The assumption that the potentials are static is approximately 
valid in a high temperature plasma of massless quarks and gluons in the following 
sense: As T — > oo, the effective coupling g — > (albeit very slowly). The perturbative 
Debye screening mass /i ~ gT limits q± ~ gT. The typical thermal energy E T ~ 3T of 
the plasma constituents is therefore large compared to /i. Consequently, the average 
energy loss per elastic collision, — q° — q z ps q\/2E T oc (? 2 T, is ~ g times smaller 
than the average transverse momentum transfer. 

Because we are interested in relatively low momentum transfer scattering (Aq C d <^ 
q± ~ gT <T), the spin of the partons can be neglected. The jet parton is allowed, 
however, to be in an arbitrary d- dimensional representation of SU(N) with genera- 
tors, T a , satisfying T a T a = C 2 \ d . 

The Born (color matrix) amplitude to scatter from an incident four momentum 
Pi-i to pf in the potential centered at Xj is then given by 

M^ft-O = 2n5(p° -plJA^e-^ , (7) 

where q« = p» — Pi-i, and is shorthand for 

A,(q,) = = -2^E T a ^(q J ) . (8) 

The differential cross section averaged over initial and summed over final colors of 
both projectile and target partons reduces to the familiar form for low transverse 
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momentum transfers: 

dai/d&iRiCi *™ , (9) 

where the color factor is 

a = ^-Tr{T a T b )Tr{T^) = C 2 C 2l /d A . (10) 

For SU(3), 2Ci gives the usual color factors 4/9, 1, 9/4 for qq, qg, gg scattering respec- 
tively. In our notation, the angular distribution is given by 

dvi/dCli = -i-Tr |^(q,)| 2 /(47r) 2 . (11) 



2.2 Sequential Multiple Scattering 

Our main assumption for computing the multiple elastic scattering amplitude is that 
the scattering centers are well separated in the sense L, = z i+ i — Zi ^> pT 1 . In a quark- 
gluon plasma at very high temperature T, l//i ~ 1/gT and the effective qg scattering 
cross section from eq.(^) is a ~ 27ra 2 //x 2 ~ g 2 /87rT 2 . Given a Stefan-Boltzmann 
density of partons p ~ 5T 3 , the mean free path is A ~ 5/(g 2 T) ^> d ~ 1/gT for 
g 1. Hence, Li ^ d is satisfied at extreme temperatures at least. 

The dominant Born amplitude for coincident sequential scattering with target 
partons from i to j without radiation is then simply 

Mjifapi-i) = J ^■■■y^M J (p J ,p J _ 1 )tA(p J _ 1 )---iA(p l )M i (p llPi _ 1 ) , (12) 

where A(p) = (p 2 — m 2 + ie) is the intermediate jet parton propagator. Amplitudes 
involving backscattering are suppressed at high energies because of the limited mo- 
mentum transfer that each potential can impart. Because of the energy delta function 
in eq.(0), the integrations over the intermediate p\ set all of them to E and lead to 
a conservation factor 

8(ji) = 27r5(p° j -pl 1 ) . (13) 

Therefore, 
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A 



k{Pk ~ Pk-lj 



P? 







n 2 

Pfc 



ie 



,(14) 



where Pq = (Eq — m 2 ) 1 / 2 « £7 , an d the product is path ordered from left to right 
with decreasing index k. Rearranging the phases in terms of the separation vectors 

R fc = L k e z + r k = x fc+i - x fc , (15) 

we can write 

Mjt^Pi-i) = SQiy-^ie^-^I^p^) , (16) 
where the reduced amplitude is 

iji(pj,Pi-i) = / \ II ^2 — „ 2 , z r ~ p^ 1 ) ' " a *(p* - p*- 1 ) • ( 17 ) 
J y k=i y Z7[ ) m) ~Pk + te j 

Because of the assumed ordering, L k > 0, and the integrals over p zk can be evaluated 
by closing the contour in the upper half plane, setting the intermediate jet legs on-shell 
with 

P* = (Po ~ Plk) 1 ' 2 « Po - plJZPo ■ (18) 

The singularities of the A k at p zk ~ Po + Kvlk + A^) 1 ^ 2 can t> e neglected because they 
leave very small residues cx exp(— fiL k ) given the assumed large separation \ih k 3> 1 
between scattering centers. Therefore, 

f | jf_1 d 2 p± k e +ipfc Rfc I 

iji(pj,Pi-i) ~ y < n f27r)2 2P f ^i(pj - Pi-0 ■ • ■ ^<(p< - p*-i) > (1 9 ) 



with p zk given by eq. (|IBD . Note that in the high energy limit Tjj survives in spite of 
the l/Po residues because A k oc Eq due to the vector nature of the coupling. Also 
the ordering of the potentials in eq. fll4]Jl9D in decreasing order of the index cannot 



be permuted in the non-abelian case because of the non-commuting color matrices in 
the Ai. 



8 



From eq.(|T9|) we can derive two interesting limits. One is the semi-classical (large 
angle) cascade limit, and the other is the eikonal (straight line) limit for multiple 
small momentum transfer scatterings. The first case illustrates how the transverse 
momentum integrations can decouple resulting in a factorized form of the multiple 
collision amplitude and is discussed in Appendix A. The second limit is however 
physically more relevant and is considered below . 

2.3 Color Algebra 

To make explicit the color algebra we write 

M ji = {a j ---a i )M;r ai , (20) 

in a shorthand notation where 

(ay-ai) =T a ^ ■■■T a * , (21) 

and we adopt the usual summation convention over repeated indices. For the color 
neutral ensemble under consideration 

m-i Mv;') (V '---(vn,,, . (22) 

Hence, 

(M;n{M-r'^) cx n jTt{t^t<) ^ n (s aka ,c 2k /d A ) . (23) 

k=i ™ k=i 

Given eq. (|23"D , the color factor associated with the jet parton is given by 

C{i, J ) = - d Tr{{a l ---a 3 ){a r --a l )) = Cl l+l , (24) 

as obtained by repeated use of the basic (aa) = C^ld relation. This is simply the 
product of the color factors, C2, occurring for each isolated collision as in eq.(|l0|). 
Therefore, even though the amplitude eq.(pO|) does not factor in color space, the en- 
semble averaged coincidence cross section does factor for large angle elastic scattering 
in a locally color neutral ensemble as shown in appendix A. 
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2.4 Eikonal Limit 

For a high energy jet, the coincidence amplitude is dominated by small angle scat- 
tering. In this case, we change variables to = Pi — Pi 1 with p zt as P ?a E 
approximately fixed, and write the total momentum transfer as 



Qji = £qj = Pi - Pi-i • ( 25 ) 

l=i 

The coincidence amplitude can then reduces to 



Mji(pj,Pi-i) « 5(i0(%---^)(-^r m 2^o / |n|^f^- x ^r(q fc )| 



x(27r) 2 5 2 (Q ±J ,-^q ±z ) . (26) 
Note that the dependence of the phase on the Zk can be factored out with 



<lzkZk ~ (Pzj - ^o)^ + (Pq - Pz.t-l)^ • (27) 

fe=i 

This phase is important only for off-shell amplitudes with p z j ^ P or p z ^-\ ^ Pq. 

To average over the transverse coordinates x fe , we employ the frozen target ap- 
proximation taking the initial and final wavefunction of target parton, k, to be 0fcj(x&) 
and by faffak) respectively. The amplitude to leave the target in a specific final state 
is obtained by replacing the phase factors by transition form factors 



e -iq fc -x fe 



F-Mk) = I d 3 x fc ^ / (x fc )e- iqfc - Xfc fci (x fc ) . (28) 



After squaring Mji, we must sum over all final states <pkf- For scattering in a chaotic 
or thermal bath we must also average over an ensemble of initial (pki- 

A simplification is possible in the high energy limit when the energy and longitu- 
dinal momentum transfers are small, and they can be neglected or replaced by their 
average values in 5(ji) and F-j. In that case, closure ( J2f 4>kf4>*kf = 1) can be applied 
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to the sum over final states. Averaging in addition over initial states with probability 
p(i), the average squared amplitude contains factors such as 

E^E^q^^-q'^H / d 3 x k e-^-^^p k ^ k ) = Th ( qxk - q ' ±k ) , 

i f J 

(29) 

where pfc(xfc) = Z)iP(^)|0fco( x fc)| 2 is the ensemble average density distribution of target 
parton k. In eq. (p9|) note that Tfc(q^) is just the Fourier transform of the Glauber 
thickness function 

T k (x± k ) = j dz p k (z,x ±k ) , (30) 

which is the probability per unit area of finding parton k at a transverse coordinate 
x_|_fc. For a broad z distribution, it would appear that we may have violated the 
assumed pL 3> 1 assumption. However, in the m! different z orderings of the centers, 
one of the previously neglected backscattering amplitudes becomes dominant and 
after relabelling the dummy indices the same result is recovered. The only essential 
assumption is that the mean free path is long compared to the range of the potential. 
With eq. (|29"D , the ensemble average of the squared amplitude is proportional to 



mM-nP^l 2 ) oc / In ^ 2 fc f 2 ^ 2 fc - oi*) V(q*) W(ql))'| 

x(2tt) 4 5 2 (q ±j , - Eq±*)<s 2 (Q±i* - XX*) • (si) 

l=i l=i 

Only diagonal color components survive because of the color neutrality condition 
eq . (|23|) . If the transverse coordinates are distributed over a radius, R ^> p^ 1 , then the 
Fourier transform of the thickness function will limit the difference, |q_i_fc — q^fcl^V-^- 
Because V k k (q k ) varies slowly on a scale 1/R <C p, we can therefore approximate 
(^4 a *(qA:))* ~ (^4 afc (qfc))* m t ne integrand. The q' k integrals result therefore in a 
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multiplicative geometrical factor 



/ III ra T *(^ - ^) 1 W V E(Hi - «^)) = Jd 2 bf[T k (b) . (32) 

^ k=i \ ' ) l=i k=i 

With the above simplification, the ensemble averaged coincidence cross section to 
scatter sequentially with partons from % to j reduces to 

da ji /(fQ±j i = J d 2 b J Yl {d 2 q± k T k {b)da k /d 2 q ±k } 5 2 (Q ±ji - J^qjj) • (33) 



k=i l=i 



This is recognized as the classical Glauber multiple collision limit, with p k (h) = 
J T k (b)da k being the probability of scattering off center fcona classical trajectory at 
impact parameter b. 

3 Induced Soft Non-Abelian Radiation 

We turn next to the inelastic amplitudes for induced radiation of a gluon with color 
c and light cone momenta and polarization 

fc" = (uj,k Z} k ± ) = [xP + ,k 2 ± /xP + ,k A _] 

^ = (eo,-e ,ei) = [0,2e± ■kJxP + ,e ± ] . (34) 

Light cone coordinates are denoted here by square brackets, [k + , k~, kj_], with /c ± = 
uj ± k z = k 2 L /k T . We chose the two physical polarization states for on shell (k 2 = 0) 
gluons to satisfy both ek = and en = with n M = [0, 2, 0j_] in terms of two 
orthonormal e±. Thus e° = e± ■ \t±/(uj + k z ). In light cone coordinates, the incident 
jet parton has p$ = [P + ,m 2 /P + ,0] with P + pa 2Eq. We focus on the soft limit 
defined by x 1 . First we consider the amplitudes for radiation from the high energy 
parton lines. Then we show that the three gluon amplitudes essentially cut off the 
soft dk±/k± spectrum at k± ~ /i. Note that the induced bremsstrahlung associated 



with a single isolated collision was derived in pQCD in Ref. ||19||. Our interest here is 



on the induced radiation pattern associated with multiple collisions. 
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3.1 Radiation from Internal Jet lines 

The amplitude to emit a gluon with color c from the jth intermediate jet line during 
sequential scattering with target partons from 1 to m is 

M%i(k,p m ,Po) = J 4^M mij+1 (p m ,pj - k) 

x {iA( Pj - k)(-2ige Pj T c )iA( Pj )} M hl (p ]}Po ) . (35) 

The integrals over p°- set their values to E because of the 5(1 j) in Mj tl . However, 
for all the subsequent intermediate lines in M m j + x the energy is shifted from E to 
Eq — u, and a new overall energy conservation factor arises 

5(ml) oc 5(p° m -E + u) . (36) 

These shifted energies change the classical momenta in subsequent legs (k > j) to 
P fc = -FLRfc with 

P u = ((E -u) 2 -m 2 )^ 2 ^P -u/v . (37) 

where vq = Po/E is the speed of the incident parton. 

To perform the p Z j integral, it is convenient to split the two propagators using 

2e Pj A( Pj - k)A( Pj ) = ^|yAfe - *) - > ( 38 ) 

which is valid for on shell radiation since k 2 = and ek = 0. The contour integral 
over p Z j can then be performed as discussed in Appendix B. 

With ( p.23|j!24| ) from Appendix B, the radiation amplitude eq.(|35|) reduces in the 
kinematic region i«l and k± <C to 

M%l « g5(ml)e- iujt ™e~ ip ™- Xm e +ip °- Xl [ f * e i( - p ^ n ^ p2 ^ 

J (27C) 2Pq 

x (]^ e+ikX3 - ]^ e+ikXj+1 ) I m , + i(Pm,V 3 )T c I hl (p 3 ,p a ) . (39) 

Note the appearance of the phases kxj and kxj + i where x^ = (tj,x.j) and the inter- 
action time, tj, is the classical transit time along the path from Xi to x, as defined in 
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eq.( 121 ). In addition, the current element in the brackets involves on-shell momenta, 
Pjyp'j, defined by 

f 3 = (E , (P 2 - (P ±j + P±j + P±i) , (40) 

f 3 = (E - to, {Pi - (P ±j + px,) 2 ) 1/2 , Pjlj + P±j) ■ (41) 

Note finally the order in which the color matrix T c appears above. 

For radiation with k± <C fi, we can factor out of the integrand a current element 
proportional to ePj/kPj. To see this, note first that for a high energy on-shell parton 
with p — [(1 — 5)P + ,m 2 ± /(l — 5)P + ,p±] and radiation with kinematics (|34|). 

kp (k ± - x P± /(l - 6)^ + x 2 m 2 /(l - 5f k\ X L ^ X ' 1 j 



This approximate independence of the current element on p^ allows us to factor out 
( f42|) from the integrals for xm± <C However, a more general expression can 
be factored out that is valid also in the high momentum transfer limit. For fixed x, 
and Eq — > oo, the momenta pj and pj are approximately fixed by geometry to be 
Pf « P n (l, R;), and therefore 



(43) 



Since this expression also reduces to ( [4"2] ) when Rj points along the jet direction, the 
final factorized form of the amplitude for intermediate line radiation becomes 

iVC'i » ^ ( e te ' W^ 1 ) (a m • ^ (44) 

In this expression, the j independent phase factor, exp(— iiot m ), was discarded. 

We note several points in connection the amplitudes for induced radiation from 
the internal jet lines given by eq. (|4"4|) . 

1 . The approximate factored expression holds both in the large angle cascade limit 
considered in Appendix A and straight line (Eikonal) limits of M m \ as long as 
i < 1 and k±_ <C //. 
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2. Just as in the QED case, these intermediate line amplitudes vanish for radi- 
ation with formation length significantly exceeding the separation of adjacent 
scattering centers. In particular, for fixed x = k + / P + ~ uj/Eq <C 1 the phase 
factors cancel in the k\ — > limit: 



oc k(xj + i — Xj) = uLj(l/vo — cos9) — kj_ • r±j « Lj/r(k) — k± ■ r±j — > , 

(45) 

where r(k) is the formation time (length) from eq.(|l]). 

3. Unlike in momentum space, there is no factorization in color space. The color 
matrix for the amplitude without radiation is modified by radiating a gluon of 
color c after the jth interaction by the insertion of a T c matrix: 

(a m ■ ■ ■ ay ■ ■ ai) — > (a m ■ ■ ■ caj ■ ■ ■ a x ) . (46) 

4. Even for i < 1 the condition leading to the approximate form in ([E3) breaks 
down for very large number of collisions since the random walk in transverse 
momentum space leads to a growing (p± m ) oc mix 2 . Thus for large m the above 
approximation is only valid in a restricted x region 

1 k ± 1 

x < — < —— . (47) 

Vm /i V m 

5. We must also add the amplitudes for radiation from the initial and final lines 
to the above amplitudes for radiation from internal lines. These external line 
amplitudes are 

AC = M ml {p m ,p -k){i^{p -k)){-2igtp T c ) 

« -9^e ik ^(a m -..a lC )M^r ai 
kP 

MZ = {tA{p m )){-2tgep m T c )M ml {p m ,p ) 

~ g^e^(ca m ---a 1 )M^ 1 - ai , (48) 

where we defined Pq = Pq, P£ = p!^ for notational convenience, and we again 
discarded the common phase factor, exp(— iut m ), as in eq. 
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The sum of all the amplitudes (|3L|S) can then be written in a suggestive form 
M^(fc,p m ,p )^^[eJ(A;)]: m ... ai M^- ai (p m ,p ) , (49) 
where the effective color current element analogous to (0) is 
M*)]L-ai =wY, eikXl [^{am-'-cavax) - j^-(a m ■ ■ ■ a lC ■ ■ ■ a{)j , (50) 



In the abelian case, the matrices in (• ■ •) are set to unity, and Eqs. (|49| , pOD reduce to 
the soft radiation formulas of QED. In the abelian case the current actually vanishes 
as a power of x in the x — > limit on account of fl4"2|) . In the non-abelian case, the 
non-commutativity of the generators leads, however, to a non-vanishing current in 
(|50D even for x = 0. 

We emphasize that the approximate form of the effective current in eq. (]50|) is valid 
only in a restricted kinematic domain (|47D. Gauge invariance requires the absence of 
induced radiation associated with collisions without momentum transfer ]19), i.e., for 
Pi = Pi-i. On the other hand, the contribution from scattering at I in eq . (|50|) is non- 
vanishing if [c,a] ^ 0. To recover full gauge invariance, of course all the amplitudes 
involving three and four gluon vertices as well must be added to the above result. 



However, since the domain of applicability of eq. (|50|) shrinks to zero as q±j — > 0, it is 
consistent with the gauge invariance requirement in the kinematic domain indicated. 

For the physically most interesting eikonal limit, the effective color current reduces 
for x^/m «C k±/ fi <C 1 to 

f ■ k m 

M^L-a, « 2z ff ^ p ^^e^(a m ...[c,^]...a 1 ) . (51) 

K ± i=i 



The lower bound on the domain of applicability comes from (|47[) . For k± < x/iy/j, 
the effective current has components in the directions p±j that cannot be factorized 
out of the elastic amplitude. We note that (|5lD can also be derived directly from 
(p5|) using the eikonal form ( p6|) for M m J+1 and M,^. The above derivation has the 
advantage that the connection to the familiar abelian case is made more transparent. 
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3.2 Radiation from Internal Gluon Lines 

In the previous section we concentrated on the amplitudes for induced radiation 
in which the gluon is radiated directly from the fast jet lines without final state 
interactions. Here we consider amplitudes involving one or more three or four gluon 
vertices in which the gluon scatters with one or more of the target partons before 
emerging with kinematics given by fl34"|). In the high energy limit, the amplitudes 
involving the four gluon vertices can be generally neglected because those vertices are 
momentum independent and a contact interaction with two widely separated target 
partons is small. The amplitudes involving three gluon vertices can be classified by 
the number of such vertices and the indices of the scattering centers to which one 
of the gluon legs is attached. The simplest of such amplitudes, denoted by G?, 
corresponds to a (possibly virtual) gluon emitted by the jet between the Zj-i and 
Zj+x with that gluon scattering off the target parton at Zj. For 1 < j < m, 
dp^ dp^ 

(2?r) 4 gfe Mw 'J +1 (P m '^'^ 

(52) 



where with qj = pj — Pj-\ the single three gluon amplitude in the Feynman gauge is 
G^Pi^-i) = (-igT b (p j . 1 + Pj ) a ){-iA(q j )) 

xi-gf^A^iq^-qj-k^^Af^ + k^ik) . (53) 
The external field at j in our case is 

A^(g)=/°27r5( g °)^(q)e- iq ^ , (54) 



and the three gluon tensor is 

A a ^ip 1 ,P2,p 3 ) = (P2 - Ps) V 7 + (ps - Pi)V a + (Pi - V2Vg aP ■ (55) 

The amplitudes with multiple three gluon vertices correspond to multiple final 
state interactions of the emitted gluon. Because the centers are assumed to be far 
apart, the intermediate gluon lines in those amplitudes are set on shell (by the cor- 
responding dp z contour integral). Those amplitudes therefore describe final state 
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cascading of the gluon in the medium. In QED these type of amplitudes are re- 
placed by higher order Compton like amplitudes. For gluons, final state cascading 
will broaden their final k± distribution and induce further gluon showering. However, 
for the problem of the energy loss of the incident jet, which is our primary interest 
here, the final transverse momentum distribution of the rescattered gluons is not im- 
portant. In addition, in the soft limit k± <C /i each triple gluon vertex gives rise to a 
factor 0(k±/fi) smaller than the corresponding internal jet line radiation amplitude, 
as we show below. We therefore concentrate here only on the amplitudes involving 
one three gluon vertex given by (|52|) . 

As these amplitudes only arise in the non-abelian case, we can simplify the deriva- 
tion by evaluating (^) in the eikonal limit. As shown in Appendix C, we find in this 
limit that 

m / j2 \ 



fc=i 



(2tt) 2 5(q ± - kL - £ an) { 2 ^' /T(fc)! 7^rk^^ m • • • t c > °A • • • ai ) 



k=l 



(56) 



We note that corrections to the effective current element in the {} brackets arise for 
large j in the region k± < xp±j ~ x[i\Jj from terms neglected in the vertex function 
in eq.( |133[ ) of Appendix C. In eq.(^j), Qx = p± m is the final transverse momentum 



of the jet parton. 

Summing these three gluon amplitudes and adding the radiation amplitudes from 
the jet lines in Eqs.([|^,[5l]), we obtain the total amplitude for m-fold coincidence 



scatterings together with soft radiation in the eikonal limit: 

M^(k,p mjPo ) « 27t6(pl-E + u J )(-zg) m (2E )jfl ^|^| e ^— ^ (q±A 

771 

x(2n) 2 S(Q ± - ki - E OiJi) {-^ ■ Jl-aJk; {q ±i })} 

k=l 
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(57) 

where the effective color current is 

m 

J^-aJk; {q ±J }) = J2e^) m(am . . . [Cj a .] . . . 0l) , (58) 
i=i 

and the elementary current elements are 

* w=2 *(S + i^&) ■ (59) 

The following points should be noted in connection with the above results: 



1. For the case m — 1, eq. floq) reduces to the result derived in [19 



2. For k± <C g^j, the three gluon amplitudes can be neglected as noted before. 
Therefore, in this limit eq.(|58|) reduces to eq.(|5"l"l). 

3. However, for very small k± < x/ixi corrections to eq.(^) arise as can be seen 
from eq.([|^). In particular, the singularities at k± = is regulated on a scale 
x/i, where \i is the (dynamic) mass of the jet parton. 

4. The singularity at k± = q^i is a non-abelian feature due to induced radiation 
along the direction of the exchanged gluon. It is regulated by the gluon polar- 



ization tensor in the medium 20 . 



5. The approximate color current is strictly valid only for k± fi and a; < 1. 
However, eq.([59|) shows the general cancellation of amplitudes for k± ^> q^i 
that limits the induced radiation from an isolated scattering to k± ~ jj,. 



6. For q±i = 0, the current element ji vanishes in accordance with gauge invariance[19 



7. The phase factor in eq.(|58D is independent of the transverse coordinates x^j in 
the eikonal limit. The transverse phase factors, exp(— zk^ -x_|_j), associated with 
each isolated collision are spread over the net elastic phase factor exp(— i J2k Q-Lfc' 

X_Lfc). 
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8. Finally, corrections to eq.(|59|) also appear in kinematic domains where either the 
radiation formation time, r(k), or the intermediate jet lifetime, 2a>/ (q_i_j — kj_) 2 , 



is on the order of the mean free path due as shown by eq. (|137|) in Appendix C 



Eqs.( |58"| , |59|) are the main result of this paper from which we derive next the non- 
abelian LPM interference effect in the eikonal limit. 

4 The Radiation Formation Factor 

The spectrum of soft induced bremsstrahlung associated with multiple scattering in 
a color neutral ensemble eq. ( f2~3|) can be computed from eq. fl58|) using steps similar to 



those leading from eq.(|2~6"D to eq . (|33|) . Analogous to eq.(|]), we find that 



- ^ TrJ2\[eJ(k,{cixi})} c ai -a m l ■ (60) 



d 3 n m 1 1 

~¥k~ = 2(2vr) 3 C^d 



Recall that C™ is the color factor for the coincidence scattering cross section without 
radiation from eq.(|24|) with C 2 and d being the second order Casimir and dimension 
of the SU(N) representation of the jet parton. As in eq.(|31|) the assumptions of 
color neutrality and that the transverse distribution of target partons is much wider 
than are essential to obtain the above diagonal form in color and q_i_j labels. 
Note that the squared current involves a sum over repeated color indices, c, a,, and 
the trace is over the resulting sum of products of color matrices that we consider in 
detail below. The average, denoted by large brackets, above is over the transverse 
momentum transfers, q_i_j, and given in our case by 

r ( m iflrftn ~1 

(/w)=/jn^j/M. (6i) 

We are mainly interested in comparing the induced spectrum for m > 1 to the 
radiation spectrum from a single isolated collision |19| : 

d 3 n! _ / a s C A qj 
U d?k ~\ vr 2 k\ (k ± - q ± ) 2 
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The average in this case is over a single scattering as obtained from eq.(|6TD by setting 
m = 1. The + label on the average indicates that the infrared singularities are regu- 
lated in a quark-gluon plasma by the dynamically generated masses of the initial jet 
and radiated gluon as discussed in the previous section. Note that this mechanism 
is different from the regulation of infrared singularities considered in Ref.|19|] due to 



form factors arising in collisions of color singlet hadrons. In that case, interference 
amplitudes associated with radiation from different hadronic constituents cancel both 
the k± = and the = q_i_ singularities. In a quark-gluon plasma at the pertur- 
bative level the quasi-particles are color non-singlet partons, and the singularities 



are regulated by medium polarization effects] 20 1. Nevertheless, in both cases similar 
expressions arise at the end. 



With eqs.( |60| , |62D , we can define the "radiation or color formation factor", C m (k) 

via 

The result can be expressed as 

i ml i— 1 \ 

= k + 2ifeE£7«eM/ T «F i ,(A) , (64) 

in terms of color coefficients, 

Qj = Tr(a m • • • [c, aj • • • a x a x ■ ■ ■ [aj, c] ■ ■ ■ a m ) , (65) 
and current correlation functions 

Fr J {k) = {j l -j ] )/m 2 ) , (66) 

with ]i given by eq. ( |59|) . 

4.1 Color Coefficients 

The color coefficients in eq. ( |55| ) can be computed by repeated use of basic SU (N) 
relations for sums of products of generators: 

aa = C 2 ld , [a, b]a = — , 
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aba = {C 2 -C A /2)b , [a, b] [b, a] = C A C 2 l d 



(67) 



The diagonal coefficients can be seen to be identical to the normalization factor in 
eq.®: 

C u = Tr(C™-%a l ]C i 2 - l [a i) c]) = C™C A d . (68) 

The off diagonal j < % coefficients are, on the other hand, 

Cjj — Tt(C™ ' [c, ctj]cij_^ • • • a,jC 2 c]ctj-|_i • • 'dj) 
= _^ic^-*+i-i(c 2 _ 9A)i-i-^c 2 d 

= , (69) 

where 

Ca__ j N 2 /(N 2 - 1) for quarks with C 2 = C F , , 

r2 ~ 2C 2 ~ \ 1/2 for gluons with C 2 = C A ' ^ 

The radiation formation factor is therefore 

C m (fe) = ra - -^fle f U(l - r^-^.^e^^ . (71) 

1 - r 2 i= l j=1 



For a single isolated scattering, of course, 

C x {k) = \ . (72) 

For multiple scattering in the Bethe-Heitler limit, corresponding to Ljj 3> t(A;), the 
phase factors average to zero, and the intensity of induced radiation is simply additive 
in the number of scatterings, i.e., 

C m (k) « m if Lij > r(k) for all i > j . (73) 

In the deep LPM limit, where r(k) ^> L^, the destructive interference pattern summa- 
rized by the formation factor depends on the form of the current correlation functions. 
It is amusing to note that that the negative sign leading to destructive interference in 
eq . (|7lD arises in QCD from the color algebra, eq. (|69|) , in contrast to QED where the 
destructive pattern in eq. (|2|) arises from the opposite sign of contributing momentum 
space amplitudes. 
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4.2 Current Correlation Function 



In order to investigate the formal structure of the current correlation function, eq. flBE| ) . 
we evaluate 

« (l + /f(*i))V*i ■ (74) 

where in terms of the transverse vector, J = (qj_ — kj_)/|q_|_ — k^| 2 

H(kl) = (krJ^kr ^^ l . (75) 

The approximate independence of the current correlations on the indices, i,j, is only 
valid in a kinematic region kj_ > xp±i recalling eqs.(f42|). In QED this restriction is 
severe because the leading term, e± ■ k±/k\, from eq.fl42|) cancels, and the photon 
spectrum is peaked at x ~ 1. In QCD, on the other hand, the radiated energy 
fraction, xdn/dxd 2 k±, is approximately independent of x from eq.(^). Therefore, 
unlike in QED the regime x <C 1 is relevant in the case of QCD. We find below that 
the induced radiation is indeed limited to x < \[i 2 /E 1, and thus the above 
approximation is justified. This approximation cannot however be extended outside 
the soft eikonal limit. For moderate x < 1 it clearly breaks down especially because 

grows approximately linearly with i due to multiple scattering. For the general 
case, the exact current element must be used and the correlation function must be 
computed from a solution of a transport equation, as first done by Migdal[|l] for QED. 
We limit the discussion here to the soft eikonal regime. 

The diagonal correlator, in the same limit is proportional to the invariant gluon 
distribution from a single collision: 



(bll 2 ) oc 



+ 



|q± • 



>i|q± 



2H{kl)+H 2 {kl))/k] 



(76) 
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which involves a second function 

HAkl) = ,i(|J| 2 ) + .( I _^- F ) + . (77) 

The approximate correlation function in the soft eikonal limit is therefore independent 
of i,j and given by 

Because (|J — (kj_ • J)/kj_\ 2 ) > 0, note that H 2 > H 2 and consequently the current 
correlation function is bounded: 

< F(k) < 1 . (79) 

The upper bound is approached for k\ <C (qj_) ~ fi 2 . In that soft region both 
H « H 2 (x k 2 ± /(q 2 ± ) < 1, and 

Fijik) « F(0) = 1 . (80) 

The lower bound is approached in the opposite limit, k\ ^> (gi), Formally H ~ 
— 1 — while if 2 ~ 1 + 3(5j_)/fcj_, and consequently in that limit 

oc («^>/fci . (81) 

The exact form interpolating between these limits depends of course on the proper 
inclusion of polarization effects in the medium. 

4.3 The Factorization Limit 

For fixed x ~ oj/Eq <C 1 and k± — > 0, Fjj « 1, and the formation length, r(k) = 
xP + jk\ becomes much longer than the separation of the scattering centers. In this 
case, the phase factors can be set to unity. With the help of 

m i ~ 1 1 — r f 1 \ 

EE(i-^r~ ™--(i-(i-^n , (82) 

i=i j=i r 2 \ r 2 / 
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the radiation formation factor reduces to 

lim C m (k) = C«=l(l-(1-T*r) 

/ 1- l/iV 2 (l-l/(l-iV 2 ) m - 1 ) for quarks 
V 2(1 - l/2 m ) for gluons 

.(83) 

Note also that for a given r 2 = Ca/2C 2 , approaches l/r 2 , independent of the 
number of collisions as m — > oo. This is the Factorization limit, in contrast to the 
additive Bethe-Heitler limit. The saturation value depends on the SU(N) represen- 
tation of the jet parton and causes the Ca factor in eq. (|62"D to be replaced by 2C 2 . 
Interestingly, for quarks in the fundamental representation the destructive interfer- 
ence is so effective that for N = 3, the final radiation intensity in the kj_ -C (q±) 
region after many collisions is even slightly less, l/r 2 = 8/9, than for a single isolated 
collision. However, for incident gluon jets, the induced radiation approaches twice 
that from a single collision. For exotic hybrid partons in very high dimensional rep- 
resentations of SU(N), the suppression effect in fact disappears altogether for fixed 
m as Ca/ttiC2 — > 0. This dependence of the LPM effect on the representation of the 
parton is a specific non-abelian effect in QCD. 

4.4 Ensemble Averaged Formation Factor 

For < k± <C fi, we can write 

C m { U , kj =C° m + —^Re E E(! - ^'(1 - F^e^y^) . (84) 

1 - r 2 i= l j=1 

To see analytically how C m interpolates between and m as a function of the ratio 
of the mean free path to the formation time, we average now over the interaction 
Z{ points according to linear kinetic theory. Because we restrict the discussion here 
to the eikonal case, the complication due to the full 3D transport evolution can be 
neglected. In linear kinetic theory the longitudinal separation between successive 
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scatterings, = z i+1 — z i: is distributed simply as 

exp(-L,/A)/A , (85) 
where A = (cr pt>o) _1 is the mean free path. Therefore, 



/ e i(*i-*i)M*)\ ~ ( \ ^ ( 86 ) 

{ ' \l-iX/r(k)J ' 1 J 

where r(k) = 2xEo/k' 2 L . 

The sum in eq . (ffl|) can then be performed since in this soft eikonal limit Fij(k) ~ 
F(k) ~ 1 In that case eq. (|82]) can be used with the replacement 



The resulting radiation formation factor reduces to 
1 + X 2 {k) r 2 \{l-tx{k)) 



2 



1 - r 2 



1 - ir 2 x(k) 



where the dimensionless function controlling the non-abelian LPM effect is 

m * = »L = . (89 ) 

T{k)r 2 2xr 2 F l-a cosh(yJ 

The last form is in terms of the rapidity, y of the radiated gluon (u = k± cosh y) . 
Note that eq.(|88D satisfies all the previous limits considered (m — 1, k± = 0, A = oo). 

For moderate large m the term proportional to (1 — r 2 ) m term can be neglected, 
and the formation factor simplifies to 

This illustrates clearly how the radiation formation factor interpolates between the 
r = and r = oo limits as a function of the dimensionless variable x- However, it 
also shows that, through the dependence on the current correlation function, F(k), 
the radiation formation factor is actually a function of two dimensionless variables, 
x(k) = X/r 2 r and k\j [i 2 . Thus, both the range, yT 1 , as well as the separation, A, of 
the interactions influences the final interference pattern. 
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For quarks, r 2 ~ 1, while for gluons r 2 = 1/2. In the collinear regime, k± ^ fi 
regime, F(k) ~ 1, and a further simplification occurs. The radiation formation factors 
for incident quarks and gluons reduce then to the following simple "pocket" formulas 



m{ 1 ~ r 2 (£;) + A 2 ' 



Am\ 2 2 

CHk) « — — + - (91} 

mK ] r 2 (A;)+4A 2 (r 2 (£;) +4A 2 ) 2 ' 1 ' 

While the above interpolation formula for quarks can be extrapolated down torn — 1, 
the gluon one only applies for large m because of the extra factor of two radiation in 
that case for large m. We emphasize again the restriction xfiy/m < fci < /i used in 
deriving the above expressions. 

5 Induced Soft Radiative Energy Loss 

An important application of the radiation formation factor is to the problem of cal- 
culation of the radiative energy loss per unit length, dE/dz, for a parton passing 
through dense, color neutral matter with a mean free path, A ^> yT 1 . We need 
only the incremental increase of the induced radiation going from m to m + 1. For 
moderately large m and kj_ < fi 

dCJdm^ X 2 (k)/(l + x 2 (k)) , (92) 

with x{k) given by eq. (|89|) . Increasing m — > m + 1, the average increase of the 
interaction length is A, and thus from eq. ( |62"| , |63"D 



at? ia d (a r Aii f Xl dxa s C A xE oX 2 (x,k ± ) 

The subscript "soft" is included to emphasize that the eikonal approximation, used 
in deriving eq.(^), restricts its applicability to the kinematic region xfj,y/m < k± < 
jj,. The limits on the fractional energy loss are xq ~ k±/Eo from kinematics, and 
x\ ~ k±/p, from the above restriction. Note that we neglect the 1/y/m dependence 
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of x\ since the Factorization limit is found to be insensitive to this cutoff, while the 
Bethe-Heitler limit must reduce in any case just the sum the radiation from isolated 
(m = 1) collisions. 

The kinematic restriction k\ < fi 2 , also follows from the destructive interference 
between three gluon vertex and jet line amplitudes. This point was already empha- 
sized in jnj but was missed in ref.[|IB[], where the domain of k± integration was allowed 
to extend up to the kinematic limit, k\ < s/4 = E m/2. That lead to the erroneous 
conclusion that dE/dz oc Eq, in violation with the factorization theorems. 

The integral over x can be performed by changing variables to Xi with the result 



x i dx xEox 2 
x x 1 + x 2 



X 2 kl{tan- 1 (\ 2 k ± )-tan- l (X 2 k ± (fi/E ))) , (94) 



where A 2 = A/2r 2 . We see that in the E 3> A 2 /i 2 limit, the second term is negligible. 
Furthermore, since A/!~l/g;^>lisa basic assumption in our multiple collision 
analysis, the first term in the brackets is approximately tt/2. 

The integral over k± is also analytic and illustrates how dE/dz interpolates be- 
tween the Factorization and Bethe-Heitler limits. We find that 

dE soft /dx « la s C 2 fi 2 {L[X 2 fj] - L[\ 2f i{fi/ E )}} , (95) 

where the interpolation function is 

i n 

(96) 



L[a] = - 

7T 



1 + -z) tan 1 (a) - - 
a z a 



Note that for a > 1, I[a] ~ 1 — 4/(7ra). For a <^ 1, on the other hand, L[a] w 4x/37r. 
Since, A 2 /i 3> 1, the first term in the brackets is always close to unity. However, the 
second term depends on the dimensionless ratio, ( = A/i 2 / 'Eq. This ratio is large in 
the additive Bethe-Heitler limit and small in the Factorization limit. 

To see how eq. (p5|) interpolates between those two limits consider first the ap- 
proximate Factorization limit. We fix A 2 /z 3> 1 and send E — > oo. In that case, the 
second term in eq.(^) can be neglected and 

dE soft /dz w l -a s C 2 ii 2 w 2vr« 2 C 2 T 2 , (97) 
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where we used fi ~ gT to estimate the force range in a quark-gluon plasma at tem- 
perature T. The result is thus sensitive to the square of the radiated transverse 
momentum as suggested in 0. Note that because we are only calculating the low 
k± < fi contribution, our derivation does not allow us to calculate logarithmic energy 
dependent factors as obtained qualitatively in |J. However, up to such logarith- 
mic factors eq.(P7p demonstrates the approximate constant behavior of the induced 
radiated energy loss in the Factorization limit. 

In the other extreme limit, we fix Eq and send A/x — > oo so that £ ^> 1. This is 
the dilute limit where the mean free path exceeds the radiation formation length. In 
this case the arguments of both terms in eq. (|95|) are large, and the small difference 
leads to 

at? ia 1 n 2 4^0 e o f 2a 8 C A \ , Qa , 
dE soft /dz ~ -a s C 2 /i — r . (98) 

Note that in this limit we recover the linear dependence of dE/dz on the incident 
energy (modulo logarithms), as in the Bethe-Heitler formula. 

It is interesting to note that in the additive regime the radiated energy loss is 
proportional to Ca, as for a single scattering via eq.(^), However, in the approximate 
Factorization limit the induced radiated energy loss is proportional to the Ci of the 
jet parton. This means in practice that gluons radiate Ca/Cf = 9/4 more gluons 
than quarks for SU(3). Recall, that the energy loss due to elastic collisions for gluons 



is also enhanced relative to quarks by the same Ca/Cf factor [^T]]. For comparison, 
the energy loss due to elastic collisions from is 

47T 

dE el /dz « — C 2 a 2 s T 2 log(£ /7ra s T) . (99) 
Therefore, the total dE/dx simply scales with C-i- This scaling differs from the quali- 



tative estimates in |9j using the single scattering bremsstrahlung cross sections of |L9 



It is interesting to note that both the elastic and radiated energy loss is proportional 
to a 2 and are comparable in magnitude up to uncertain logrithmic factors. 

We emphasize that eq. (|97| ) for the radiated energy loss is only an order of magni- 
tude estimate because we have not calculated the contribution from the non-factoring 
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k± > n domain. In order to improve the estimate, all the multiple three gluon ver- 
tex amplitudes added to those computed here. That moderate high k± regime also 
requires a more careful treatment of the current correlation function as well as of the 
polarization mechanisms that regulate infrared singularities. A proper treatment of 
the above problems remains an open theoretical challenge. 

Finally, we comment on the comparison of eq. fl97|) to the bound on dE/dx derived 
by Brodsky and Hoyer [flO|| based on the uncertainty principle. For a radiated gluon 
carrying away a fraction x of the incident energy, Eq, with a given k±, the uncertainty 
in its formation length is r(k) from eq.(|I]). The induced energy loss for radiating one 
gluon is therefore bounded by 

dE/dx < (xE/r{k)) = {k]_)/2 ~ fx 2 /2 . (100) 

Our estimate satisfies this bound because a s < 1 was assumed throughout our per- 
turbative analysis. In fact, we may interpret C<i®. s roughly the probability of ra- 
diating one gluon with r(k) < A between multiple collisions. That gluon is radi- 
ated in a cylindrical phase space with approximate uniform rapidity density between 
< y < log(/iA/r 2 ) and limited k± < ji. 

6 Summary and Discussion 

In this paper, we initiated a study of multiple collision theory in pQCD concentrating 
on the eikonal limit. We calculated elastic and inelastic multiple collision amplitudes 
for a high energy parton propagating though a "plasma" of static target partons. 
We showed that the assumption of color neutrality was vital to recover the classical 
parton cascade picture in both large and small angle scattering. Our main simplifying 
assumption was that the mean free path (A ~ l/g 2 T) in the target was large compared 
to the range of the interactions (/i -1 ~ 1/gT). In particular, we showed how the 
classical Glauber scattering cross section eq. (|33|) emerges after ensemble averaging. 

The main focus of the paper was to derive eq.([57]), which shows how the sum of 
the induced gluon radiation amplitudes can be expressed as a convolution of elastic 
multiple scattering amplitudes and an effective color current, eq. fl58|) . This result was 
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derived for the soft radiation (i <C 1 and k± <C //) regime. The limitations of the 
approximations leading to that result were also carefully analyzed. We showed that 
the triple gluon diagrams can be neglected in the soft limit but are important to cut 
off the k± distributions on the scale \i. Only when k± —>■ can the effective current 
be pulled out of the multiple collision integral, and the radiation amplitude factored 
in momentum space as in eq. Q39|) . There is never a factorization of amplitudes in 
color space. However, the color neutrality condition, eq. , greatly simplifies the 
ensemble average of the squared amplitude. In addition, the off diagonal contributions 
in momentum space from J^.,,^^; { c LLi})^ai-a m (^; {^li}) drop out, as shown below 
eq . (|3l|) , if the transverse width of the target is large compared to the interaction 
range. Under these conditions, it was possible to calculate the induced radiation 
spectrum from eq. (|60|) ■ 

We defined the radiation formation factor, eq. (|64j) , as the ratio of the induced 
radiation spectrum to the spectrum from an isolated collision. That factor measures 
the suppression of induced radiation with formation length, r(k) > A, and reveals 
the non-abelian analog of the LPM effect. The novel role of the color algebra that 
leads to this destructive interference pattern in QCD was shown in Eqs. ( ]69|J7T| ). We 
showed how this factor interpolates between the saturated Factorization limit and the 
additive Bethe-Heitler limit. A compact "pocket" formula for the formation factor 
was derived in eq.(|9l|) illustrating the essential features of that interpolation. 

Finally, we applied the formation factor to estimate the contribution of soft in- 
duced gluon radiation to the energy loss per unit length, eq.(95). The result in the 
Factorization limit, eq.(^), was shown to be consistent with the uncertainty prin- 



ciple bound of [10] with a numerical coefficient, C^ctg, that had a simple physical 
interpretation as the number of induced gluons radiated in the limited phase space 
with rapidity between zero and log(/iA/r 2 ) and with kj_ < ft. Up to un-calculated 
logarithmic factors the radiative energy loss was found to be comparable to the elas- 
tic energy lossf^l]]. We also showed how the linear energy dependence of dE/dz is 



recovered in the opposite (Bethe-Heitler) limit when X^> E/^j 2 . 

Naturally, many problems need further study. Especially important will be to ex- 
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tend the derivation of radiative cross sections to the moderate k± > fj, regime to cover 
the non-factorizable "semi-hard" regime for induced radiation. We concentrated on 
the soft regime here to simplify our task. The "semi-hard" regime is however also 
complicated by the necessity of having to consider in detail the polarization effects 
that regulate = singularities and also the necessity of computing the current 
correlation functions discussed in section 4.2. The open theoretical question in this 
connection is to what extent, if any can a classical parton cascade transport model 
be constructed that correctly simulates the many subtle interference phenomena of 
pQCD in the multiple collision domain. 
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Appendix A: Large Angle Elastic Cascade Limit 

We consider here the elastic coincidence scattering in the special idealized case that 
the coordinates Xj are fixed and the energy is large to see how the classical billiard 
ball formula emerges. We rewrite the phases in eq.flTSp as 

Pfc • R/c ~ Pfc ■ R-fc — 77^-(P_Lfc — P_Lfc) 2 , (101) 

ZJr'o 

with Pfc = P R, k /R k and 

P±k = P r ±k /R k ^P r ±k /L k . (102) 

Note that p_L fc ■ r_L fc = 2p±k ■ P± k {L k /2P ) and P L k « P fc ■ R k - P 2 ±k (L k /2P ). 
Substituting eq.(|101|) into eq.flT9|) and shifting the p± k integration, we find that 



X 



i-1 

II^(Qi + Pii-Pi(i-i))l , (103) 

l=i 



where the intermediate classical momentum transfers, denoted by 

Q k = P k - P fe _! for i < k < j , (104) 

have dominantly transverse components. The endpoint momentum transfers are given 
Qj = pj — Pj-i and Qj = Pj — p^-i, and the external phase is 



A simplification occurs in the high energy when the Xj are fixed because the Q±fc 
increase linearly with E , while the transverse momentum integrals are limited by 
the oscillating phase factors to p± k <(P /L k ) 1 ^ 2 . Hence Q± k 3> p± k for high enough 
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energies, and we can expand the potentials around Q± k . This is equivalent to the 
stationary phase approximation. The region of validity of that approximation can be 
clarified by writing 



Q±k = Po(r ±k /L k - r ±fc _i/L fc _i) = P Ar ±k /L k . (106) 

The condition that Q± k 3> (Po/L k y/ 2 is thus equivalent to requiring that the trans- 
verse momentum transfer be large enough to resolve the transverse separation of the 
scattering centers: 

Q± fc -Ar ±fc >^ . (107) 

For fixed rj_ k this condition is always satisfied for sufficiently large energies, and thus 
the stationary phase integrals can then be evaluated using 

In this limit Mji reduces to the simple factorized form 



j'-i 

MAl',-1'' ~ 5(j0e + ^- x ^ J (Q,)n{e m/2 e +lQ - x ^ fe (Q fc )/(47ri? fc )} .(109) 

k=i 

Note again that in the non-Abelian case the matrix ordering from j to % is essential. 

After squaring and integrating over the magnitude of the final momentum, aver- 
aging over initial colors and summing over final, the above factorized form leads via 
eq.flllD to the classical billiard ball formula 



daji/dn = dai/dtti I J] da k /R 2 k dtt k \ . (110) 

[k=i+l J 

We emphasize that color neutrality of the medium and large transverse momentum 
transfers are essential to recover this simple cascade limit in which the direction of 
all intermediate momenta are fixed by geometry. 
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Appendix B: Derivation of Eq.(39) 

Technical details of the derivation of the radiation amplitude from intermediate jet 
lines are given here. The integration over p z j in eq.fl3"5|) can be expressed as a sum of 
two terms using eq.(^). The term containing the A(pj — k) propagator leads to 

d 3 pj < H ( ' k x 



M+ = tg5(ml)e~ ipm -* m e +ip °-* 1 J 



(2tt) 3 PI - ( Pj - kf + % e 



X ^ _ k ^ Im,j+l(Pm,Pj ~ k)TJj tl (p j} Po) . (Ill) 

The term containing the A(pj) propagator leads to 

CJ v ' J (2vr) 3 P 2 -pf + «e 

€.7) ' 

Xt— I mtj+ i(p m ,pj - fyTJj^ip^po) . (112) 

Kpj 

In the contour integral over p Z j sets pj — k on shell with 

p zj ~ fc 2 + P w - (P±j ~ k±) 2 /2^ 

« P -(p± J -k ± ) 2 /2P + (A; z -^ ) , (113) 

while in M~j it sets on shell with p Z j m P — p^/2P - In both cases the residue 
of the propagator can be well approximated by 1/2Pq in the high energy limit when 
x < 1. However, it is essential to keep track of the difference, P — P w = lo/vq in 
computing the phases. The phase in the integrand of is given by 

Pi R i - k x i+i = (Pi - k)Rj - k ■ Xj 

« Pi ■ Ri " ^r(P±i - k± - P±,) 2 - uRj/vo - k • x, (114) 

Note that the replacement uoLj/vj — > uRj/vj above is valid either when rj_j ~ <C Lj 
or P — > oo. The phase in is, on the other hand, 

Pj R-j " k *i+i « P, ■ R, - ^r(P±i - P±i) 2 - k • x i+1 . (115) 
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Shifting the p±j integration to p±j + kj_ + P±j and p±j + P±j in M^j respectively, 

M+ « gS(ml)e- iPm Xm e +ip °- Xl f ^1—^^^^% 
CJ y K ' J 27T 2 2Pn 



x \J^ e -^/"i+^)j J mJ+1 (p mj p i )r < J iil (p ilPb ) , (116) 



X 



^e- ik ^+M/ roJ+1 (p ro ,p i -A:)T c /,- 1 (p i ,po) • (117) 



Note that in M~j, the on-shell pj and pj are given by Eqs. (|40|j4l|) . Therefore both the 
off-shell amplitude I m ,j+i(pm,Pj — k) and the on-shell amplitude I m j + i(p m ,pj) in the 
integrands above are evaluated with the same shifted incident energy, E — u, and 
approximately the same shifted incident longitudinal momentum, P^ w Pq — uj/vq. 

The energy shift of I m ,j+i above leads in the high energy limit to a phase shift of 
those amplitudes relative to the case without radiation. To see this, note from eq.([jT 
that 



I m , j+ l(Prn,Pj-k) = J < J] (27l)3 p 2 _ 2 T J \ ^(Pm-Ptn-l) ' ' ' ( Pj+1 - Pj +k) 



k=j+l 



(2tt)3 PI - p\ + ie 



118) 



The contour integrals over the p zk then fix the phases to be 

p k R fc w P w L fc - p 2 Lk (L k /2P L ,) + p ±k ■ r ±k 

« P fe • R fc - -^-(p± fe - P± fc ) 2 - • (H9) 

Therefore, there is an additional phase shift uL k /vo « uR k /vQ for each intermediate 
line. This phase shift has a simple physical interpretation. Noting that the classical 
transit time between centers at x& and x& + i is 

At k ^R k /v , (120) 
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the additional phase shift is due to the time delay of a wavefront propagating with a 

frequency reduced by uj. Defining the classical interaction time at center, k, as 

fe-i 

t fc = ^At, , (121) 
i=i 



with t\ = 0, the net extra phase can be factored out from eq.( |118| ) as 



m_1 rl 2 n P ~ i 2ifc( p±fc ~ P - L '=) 



k=j+l 



(2tt) 2 2P 



x A m (p m - p m _i) • • • A j+1 (p j+1 - pj + k) , (122) 

where <p m ,j+i is the phase without radiation given via eq.( |105[ ). In eq.( [L22| ) we again 
used the condition x <C 1 in replacing the residue 1/P U by 1/Pq. 

Note that the longitudinal momentum transfer in all the potential is still small 
since all the intermediate longitudinal momenta are shifted by approximately k z . If 
in addition to i < 1, we consider radiation with k± <C //, then the arguments of all 
the potentials inside the integrand can be approximated by those in eq.fllDp. Only 
the extra u dependent phase must be kept. Therefore, in this soft limit 

I m j + i(p m ,Pj - k) « e-^-^Imj+^Pj) , (123) 

where the right hand side is to be evaluated ignoring the soft radiation via eq.(|l9|). 
Similarly, it follows that for soft radiation 

Wl(Pm>^') ~ e " M<m "' J+l) Wl(Pm,Pi) , (124) 



involving the same phase shift as in eq. (|123|) . Combining these results we obtain 
eq.(39). 

Appendix C: Derivation of Eq.(56) 

The details of the derivation of the radiation amplitudes involving one three gluon 
vertex are given here. From eq.(p6|) 

Mj-i^(pj-i,p ) = 27r < J( P ;_ 1 - i? )(a,_i • • • a 1 )(-i g y- 1 (2E )e- i ^U-^ Po >^ 
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k=l 



(27T) 



2 e" iq — ^(q ±fc ) (2 7 r) 2 5(p ±0 -_ 1) - £ q ±fe ) 



fc=i 



(125) 



M mJ+1 (p m , Pi ) = 2vr5(^-p°)(a m ---a J+1 )(-^) m - J (2£;o) f 



771 / ?2 \ 7YI 

x ' II brrl e " q±fc ' xxfc ^ afc (q± fc ) (27r) 2 5( P±m - P± , - £ q±fc ) 



4 + i V ( 2 ^ 



k=j+l 



(126) 



We have utilized above that p z0 = P Q ^ E Q and that the energy of internal lines k < j 
is E , while for k > j it is E — u> with p z k = P^ = Pq — uj/vq on those lines. Note 
that the amplitude M m J+1 differs from the case of no radiation eq. (|26|) by only a z 
dependent phase obtained by replacing P with P u in eq . (|2~7D . In the eikonal limit 
eq. (|53|) can be expressed as 



d 2 q±j 
(2tt) 2 



, * ^;"(q..,)(2^^(q., - (p ±j - p ±( ,._ 1} + kj) 



X 



-ir(k,pj,Pj- 



u 2 - (p zj - p z(i -i)) 2 - (q±j - kj 



where we used —iT b f ba . c = [c, ctj] and note that the vertex function is 



r(k,p j ,py 



AE ( Pj - Pj-i)e(k) + Aoj Pj e(k) 

-(Pj+Pj-i)(Pj -Pj-i + 2k)e°(k) 



(127) 



128) 



Inserting these expressions into eq.(^), the integrals over the variables 
give rise again to an overall 5(ml) factor eq. (|36|) and set p®_ x = E and p° = E u = 
E — to. Then as in eq. ( |1 1 1|) we integrate the contour over p Z j and p z (j-i) keeping 
only the residues at the poles of the propagators since Lj 3> d. Because of the three 
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propagators involved, A(pj)A(j)j — p 3 -_ 1 )A(p J -_i), those integrals give rise to three 
contributions corresponding to forward scattering with two of the three internal lines 
set on shell. A fourth contribution corresponding to backscattering of the gluon near 
Zj is suppressed in the kinematic range E 3> u> 3> max(k±,q±j) of interest here. 
The kinematic variables of the three forward amplitudes after the contour integration 
are then as follows: In case 1, p 2 = p 2 _ x = 0, (pj — Pj-i) 2 ~ — (q_i_j — kj_) 2 , with 
P±j - P-L(j-i) = Qjlj - k_L, and 

Pj-i = ( E o,E -p 2 L{j _ 1) /2E ,p ±(j _ 1) ) 
Pj = {Eo-UtEo-u-pljftEo,?^) . (129) 

In case 2, (pj — Pj-i) 2 = p 2 _ x = 0, p 2 « — (q±j — k±) 2 /x, and 

Pj-i = (E Q , E - p 2 ± ( j _ 1 - ) /2E , p±(j-i)) 
pj = (E - lu, E - lu + (q_Lj - k 1 _) 2 /2u, p ±j ) . (130) 

Finally, in case 3 (pj — Pj-i) 2 = p 2 = 0, p 2 _ 1 ~ (q±j — k±) 2 /x, and 

pj-i = (E , E - (q ±j - k±) 2 /2u, P_l(j-i)) 
pj = (E -u,E -LU-p 2 Lj /2E ,p ±j ) . (131) 

We assume that lu 2 ^> (q_Lj — kj_) 2 . It is remarkable that the residue of the product 
of propagators is approximately same up to a sign in all three case with 

{-l)Res{A(pj)A(jpj - p^OAfo-O) « ±- \ , (132) 

with + for cases 1 and 2 and — for case 3. Also the vertex factor turns out to be 
approximately the same in all three cases 

r«-4£ ei-(qu--k ± ) . (133) 



For a jet parton with mass, /x, the singularity in eq. ( |132|) at q±j = k± is automatically 
regulated as in eq.(^) by 

|q±i - k ± r 2 -> (|q ±j - kj 2 + xV)" 1 . (134) 
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In higher order, the inclusion of the gluon proper self energy tensor in the medium 
would replace this infrared regulator by the effective dynamic mass of the gluons. 
Only the z dependent phase, 

<\>z = -{Pu>- Pzj)z j+ i - (p zj - Pz(j-i) + h)zj - (p z (j-i) - Po)zj-i , (135) 
is found to be case dependent. Evaluating with eq.( |129| . |T30|J131|) , we find that 



!0 case 1 

Lj/rj case 2 

—Lj_i/Tj case 3 

(136) 

where Tj = 2o;/(q_Lj — kj_) 2 is the lifetime of the virtual jet line in cases 2 and 3, 
and Lj = Zj + i — Zj is the longitudinal distance between adjacent scattering centers. 
Together with the relative signs in eq.( |132j ), the phase factors in the three cases sum 
to 

We therefore find a new interference form factor, fj , that involves the separation 
distances between adjacent centers and the lifetime of the virtual jet state, Tj. For 
well separated centers, in the sense that the mean free path A = ((^+1 — Zi)) ^> Tj, 
the extra phases in fj average to zero and fj ~ 1. Also in the extreme opposite 
limit, A Tj, extra terms tend to cancel again leading to fj ~ 1. In particular, for 
= qij, fj = l. As a rough form illustrating these limits, 

fj(k) ~ 1 + 2isin(A(q ±j - k ± ) 2 /2cu) . (138) 

In the general case, fj ^ 1 reflects the effects of final state cascading of the emitted 
gluon. Another important limit is qj = or kj_ ^> qj_j. In this limit, Tj ~ T(k) = 
1/(lu — k z ), and fj reduces to 

fj w 1 + e ^ /r(fc) - e iL ^ l/r{k) . (139) 
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For radiation with formation time much less or much longer than the mean free path, 
fj ~ 1. Therefore, except in the restricted kinematic domain where Tj or r(k) is on 
the order of the mean free path, this extra interference effect can be neglected, and 
fj can be set to unity. Combining these results, we obtain eq.(56). 
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